AD- A 187  151 


,«■  tl  J|  .»«  .«  Mm*  l.«  j.tj.f  ».t  |4t| 


SECURITY  CLASSIFICATION  OF  THIS  PACE  flMJ.il  Dmim  Knlmrmd) 


REPORT  DOCUMENTATION  PAGE 


7  -  10  81 


liiiM 

warn 


m 


2  GOVT  ACCESSION  *0-1  J.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (md  Sabllll,) 


Characterization  of  Nonhomogeneous 
Poisson  Processes  Via  Moment  Conditions 


5  Type  OF  REPORt*  PERIOD  COVERED 

Technical  -  August  1986 


7.  AuTHORfJ) 


Zhao ben  Fang 


»  performing  ORGANIZATION  name  and  address 

Department  of  Mathematics  and  Statistics 
University  of  Pittsburgh 


ii.  controlling  office  name  and  aodress 

VK\  410 

“DC  4  V 


»4.  MONITORING  AGENCY  name  ft  ADDRESSES  dltlmrmnt  from  Controlling  Olllcm) 


6.  PERFORMING  org.  report  number 


nrair.i 


i.  contract  or  grant  numbers 

AFOSR-84-0113 
ONR  N00014-84-084 


CZ-  O  lA 


I).  number  of  pages 
19 


IS.  SECURITY  CLASS,  (ol  ihl,  r.pori; 


Unclassified 


IS*.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


<6  DISTRIBUTION  STATEMENT  (ol  Ihlo  Rtpoel) 


i  \i  :  >  v 


OlbTRiauTiQN  bT  AT  EMEN  T  lot  l/i*  mb»trmct  ontwfd  In  block  30,  It  dUloronl  from  Hoport) 


»®  supplementary  notes 


<*E 


'9  KEY  W0R0S  (Contlnum  on  rmvroo  aid*  It  nacoaaory  and  Idmntlty  by  block  num bmr) 

Nonhomogeneous  Poisson  Process,  independent  incremants,  characterization 
probelm,  joint  cumulant,  exchangeable. 


20  ABSTRACT  (Continue  an  tovottm  ,1dm  II  n,c,„»ry  "id  Identity  by  block  number; 

The  property  of  independent  increments  is  one  of  the  most  important  for 
defining  both  the  homogeneous  and  nonhomogeneous  Poisson  process. 

In  this  paper  we  give  two  ways  to  relax  this  requirement  and  characterize 
the  nonhomogeneous  Poisson  process  by  some  moment  conditions. 

On  result  is  that  a  counting  process  {N(t),t  >  0}  with  finite  moments  of 
all  orders  is  a  nonhomogeneous  Poisson  process  with  mean  functions 


m(t)  *  EN(t)  if  and  only  if  for  any  t^,  i  ■ 

cum  (N(t,  ) , . .  .  ,N(t.  ))  ■  min  EN(t.) 

•  *  Ki^k  1 

where  cum  (•)  is  the  joint  multivariate  cumulant. 

A  second  result  is  that  if  increments  on  any 
interval  are  Poisson  distributed  and  an  exchangeable 
condition  is  assumed  then  the  process  is  nonhomogeneous 
Poisson.  This  extends  Renyi's  (1967)  result. 


i  ^iv>n.L»i  r 


AFOSR.TR*  87  -  1  0  81 

CHARACTERIZATION  OF  NONHOMOGENEOUS 


POISSON  PROCESSES  VIA  MOMENT  CONDITIONS 


Zhaoben  Fang  1,2 


SERIES  IN  RELIABILITY  AND  STATISTICS 


Department  of  Mathematics  and  Statistics 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 


Technical  Report  No.  86-06 
August  1986 


ter) 


Accession  For 

«XS  GRAtl 
DTIC  TAB 
Unannounced 
Justification 


I 


By-— - 

Dlstrlbut Ion/ 

Availability  Codes 
Avail  and/or 
Dlst  Special 


Department  of  Mathematics  and  Statistics,  University  of  Pittsburgh,  Pittsburgh,  PA  15260. 
Supported  by  Grant  No.  AFOSR-84-Ol  13  and  ONR  Contract  N00014-84-084. 


CHARACTERIZATION  OF  NONHOMOGENEOUS 


POISSON  PROCESSES  VIA  MOMENT  CONDITONS 

by 

Zhaoben  Fang  1,2 


Department  of  Mathematics  and  Statistics 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 


Abstract 

The  property  of  independent  increments  is  one  of  the  most  important  for  defining  both  the  homo¬ 
geneous  and  nonhomogeneous  Poisson  process. 

In  this  paper  we  give  two  ways  to  relax  this  requirement  and  characterize  the  nonhomogeneous 
Poisson  process  by  some  moment  conditons. 

One  result  is  that  a  counting  process  {N (t),t  >  0}  with  finite  moments  of  all  orders  is  a  nonhomo¬ 
geneous  Poisson  process  with  mean  functions  m(t)  =  EN(t)  if  and  only  if  for  any  tj,  i  = 


cum  (N(ti),...,N(tk))  =  min  EN(t) 
l<i<k 


where  cum  (•)  is  the  joint  multivariate  cumulant. 

A  second  result  is  that  if  increments  on  any  interval  are  Poisson  distributed  and  an  exchangeable 
condition  is  assumed  then  the  process  is  nonhomogeneous  Poisson.  This  extends  Renyi’s  (1967) 
result. 
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joint  cumulant,  exchangeable. 
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1.  Introduction 

\v 


Poisson  processes  play  an  important  role  in  many  fields.  The  Poisson  process  is  one  of  the  sim¬ 
plest  counting  processes  and  is  a  building  block  for  many  other  processes,  especially  for  general 
independent  increment  processes. 


Many  definitions  have  been  given  for  a  counting  processes  to  be  a  Poisson  process  and  many 
papers  have  appeared  dealing  with  characterizations  for  the  Poisson  distribution  and  the  Poisson  pro¬ 


cess. 


Among  the  qualitative  conditions  defining  a  Poisson  process,  independent  increments  is  one  of  the 
most  important  conditions.  In  this  paper  we  attempt  to  use  other  conditions  in  place  of  independent 
increments.  This  provides  a  somewhat  different  viewpoint  for  examining  Poisson  processes.  In  addi¬ 
tion,  new  characterizations  for  the  nonhomogeneous  Poisson  process  via  moment  conditions  are 
obtained  which  might  be  easier  to  utilize  in  practice.  " 


For  clarity  we  give  a  standard  definition  first. 


Definition  1.  A  counting  process  (N(t),  t  >  0}  is  said  to  be  a  nonhomogeneous  Poisson  process  with 
parameter  function  m(t)  if 


(i)  P{N(0)=0}  =  1, 


(ii)  the  process  has  independent  increments, 


(iii)  for  0  <  s  <  t,  N(t)-N(s)  is  Poisson  distributed 
with  mean  m(t)-m(s),  i.e. 


P{N(t)-N(s)=k}  =  {-m(-)  exp  {-(m(t)-m(s))} 

k  * 


(i) 


for  k-0,1,2,...  where  m(t)  is  a  nonnegative  nondecreasing  function  of  t. 


V 
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For  the  purpose  of  exposition  we  assume  m(t)  is  continuous.  Renyi  (1967)  made  an  essential 
improvement  in  relaxing  the  condition  of  independent  increments.  This  is  given  in  the  Theorem  1 
below. 

Let  R+  =  {t,  t  ^  0}.  B  be  the  Borel  set  of  R+,  {finite  union  of  disjoint  finite  intervals  (a,bj 
on  R+},  and  I=(s,t]  be  an  interval  on  R+.  Let  ^(E)  be  an  additive  stochastic  set  function  defined  for 
Eej?’ which  counts  the  number  of  events  of  a  counting  process  N(t)  falling  in  E.  In  other  words,  £(E) 
represents  the  total  increments  of  the  process  N(t)  on  intervals  in  E.  It  follows  that  ^(1)  =  N(t)-N(s) 
is  the  increment  on  (s,t].  Let  X(E)  be  a  measure  with  no  atoms  in  E  £y*\  such  that  X(E)  <  °°  for 
E£y£  and  X(I)  =  m(t)-m(s)  where  I  =  (s.t].  Here  m(t)  is  a  nonnegative  nondecreasing  function  of  t 
which  is  assumed  to  be  continuous. 

Theorem  1.  Let  (N(t),  t  t  0}  be  a  counting  process  with  P{N(0)=0}=1.  If  for  Eej^  ^(E)  has  a 
Poisson  distribution  with  mean  value  X(E),  then  N(t)  is  a  Poisson  process  with  mean  m(t). 

Proof.  Renyi  (1967). 

If  it  is  assumed  only  that  ^(1)  has  a  Poisson  distribution  for  any  interval  I  in  R+,  Shepp  (see  Gold¬ 
man,  1968)  has  constructed  a  counter-example  to  show  that  the  process  need  not  be  Poisson.  Further¬ 
more,  Oakes  (1972)  has  provided  a  construction  of  a  counting  process  such  that  even  if  for  fixed  k  the 
counts  ^(Ij)  in  any  k  contiguous  interval  Ij  =  (tj_j,tj]  are  independently  Poisson  distributed  with  mean 
tj-tj_j.  N(t)  need  not  be  a  Poisson  process.  This  process  is  called  a  k-fold  quasi-Poisson  process. 

An  interesting  question  is  to  determine  additional  conditions  needed  to  ensure  that  a  process  is 
nonhomogeneous  Poisson  with  mean  m(t)  besides  the  condition  that  ^(1)  is  distributed  as  Poisson  with 
mean  X(I).  In  Section  2,  we  introduce  an  exchangeable  condition  which  was  suggested  by  Dr.  Z.  D 
Bai  (personal  communication)  to  replace  the  independence  condition. 

It  is  well  known  that  independence  implies  uncorrelatedness,  but  the  converse  is  not  true  in  gen¬ 
eral.  In  Section  3,  we  use  the  idea  of  joint  cumulants  developed  in  Block  and  Fang  (1985)  to  give  a 
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characterization  for  nonhomogeneous  Poisson  processes  via  moment  conditions  under  which  uncorrelat- 
edness  between  increments  implies  independent  increments.  It  naturally  follows  that  if  ^(1)  is  Poisson 
distributed  for  any  interval  then  the  counting  process  is  nonhomogeneous  Poisson. 

Both  arguments  above  relax  the  common  condition  of  independent  increments  for  a  process  to  be 
Poisson. 


2.  An  exchangeability  condition 

In  this  section  we  prove  that  a  counting  process  is  Poisson  under  an  exchangeability  condition. 
Here  we  use  exhangeability  with  respect  to  a  finite  number  of  rv’s. 

Consider  a  counting  process  N(t)  on  R+.  Let  I  be  the  interval  (a,b],  1 1 1  =  b-a  (the  Lebesque 
measure  for  interval  I).  A(I)={N(b)-N(a)=0}  (the  event  that  there  is  no  occurrence  of  the  process  in  I), 
and  let  lA(ij  be  the  indicator  function  of  A(I). 

First  we  prove  a  lemma  to  show  how  the  exchangeable  condition  works,  then  by  using  this  lemma 
we  c  btain  a  new  set  of  conditions  for  a  Poisson  process. 

Lemma  1.  If  for  any  interval  I  =  (a,b]  on  R+ 

P{N(b)-N(a)  =0}  =  e~X|11  (2) 

and  for  any  contiguous  intervals  I4  =  (tj—i.tj]  i=l,...,k  such  that  |Ij|  =  5,  8  >  0,  the  rv’s 

1  A(ij)>  i=l,...k  are  exchangeable,  then  1  a(I J  are  independent. 

M 

Proof.  From  A(Ij)  =  {N(t^)-N(ti_i)  =  0}  we  know  that  if  J  =  Iim,  im£  {l,...,k}  m  =  1,...,M 

m=l 

M  M 

then  A(J)  =  ^  A(I,  )  and  1  a(J)  =  I”I  1  A(I  )•  is  easy  ,0  see  that  ^  A(I,)  are  identically  distri- 
m=l  m=l 


buted  Bernoulli  rv’s  with  parameter  p,  where  p  =  e 

Similarly,  since  the  1  A(I,)  S  are  exchangeable  we  have 


PUa(J)=  1}  =  P{1a,i.  )=l.m=l . M} 
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=  P{lA(li)=l,  i  =  l,...,M} 

=  P{N(tM)-N(t0)  =  0} 

_  e->-l(totM]l  _  e-X(tM-to)  _  e-X8M  _  pM 

i.e.,  lA(j)  is  Bernoulli  with  parameter  pM. 

Now,  for  fj£  {0,1},  i  =  Let  B0  =  {i|fj  =  0},  Bt  =  (i  1  fj  =  1}.  B0  Bj  =  {l,...,k}. 

Card  B0  =  M0,  Card  Bj  =  M!  and  M0  +  M!  =  k.  Denote  the  complement  of  A(I)  by  AC(I).  By 
exchangeability  we  have 


P  -  fj,  i  ~  1*' ->k} 


-  p  {!a(D  ~  O.ieBo;  1  A(ij)  —  1,  ieB!> 


—  P  { ^ A(ij  -  0,  i  -  1,...M0;  lA(ij  =  1,  i  -  M0  +  l,...,k} 


=  P  {  lA‘(y  =  1,  i  -  1,...M0,  lA(ij  =  1,  i  -  M0  +  l,...,k} 


=  E 


M„ 

n  U'oj 


n  i 


i  =  1 


i  =  M0  +  1 


A(U 


E  I  1  Mq  1  k 


n  ac(i,)  n  ac.) 

i  =  1  1  =  Mo  *  1 


E  J(1  -  lMo  )  1  « 

1  u  A(IJ  o  ACi) 

1  :  _  i  1  '  *  l 


E  , 

*  •* 

1  > 

E  • 

[Mo 

I  c-ir1  Sj 

n  1  A(u } 

n  A(ij 

j  =  i 

i  =  Mo  +  1  J 

_»I  ,||  k^  .*»  .l|  .•».«»  .»>  .U  -»l  .il 


where  S}  =  £  lA(Ij }  n„.n  A(I|).  Since  the  lA(Iy 

1  <  i[  <■■•<  ij  S  M0  1  1 


s  are  exchangeable 


E  iSj  FI  1A(Ii) 

i  =  M0  +  1 


x  E 1  *  ... ,  n  Udi) 

1  <  i,  <•••<  LSM0  n  A(  V  i  =  M0  +  1 

m  =  1 


X  E  Uajpi-nAOiinAOM^ii-nAOk) 

1  <  i,  <•■•<  ij  S  M0  1 


'  ■> 

=  2  p^U  pipM 

1  <  i!  <  --<ij  <  M0  ^ 


It  follows  from  (3)  and  (4)  that 


pk*M„  _  £  (_,y-,  Mo  pi  pk-Mo 

j  =  l  J 


=  PK-M°  (1-pf 


-  n p<  Uaj  -  fi| 


We  now  state  the  main  result  of  this  section.  Here  the  usual  condition  of  independent  increments 
is  replaced  by  an  exchangeability  condition. 

Theorem  2.  Let  (N(t)  t  >  0}  be  a  counting  process.  If  the  conditions  i)  -  iii)  below  hold,  the 
{N(t),  t  >  0}  is  a  homogeneous  Poisson  process  with  mean  function  E{N(t)}  =  Xt. 

i)  For  any  interval  I  =  (a,bj  P{N(b)-N(a)  =  0}  =  e-*-111. 


VVV.V.V.V.V.  /.  '•»  -  ./-.r.  W 
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ii)  P{N(b)  -  N(a)  >  2}  =  o  (X 1 1 1 )  as  |I|  — >  0  uniformly  in  I 

iii)  For  any  finite  number  of  contiguous  intervals  Ij  =  (t;_j,  tj]  i  =  l,...k  such  that  1 1;  |  =5,  the  rv’s 

1  A(ij  i  =  1 . k  are  exchangeable. 

Proof.  By  the  proof  of  Theorem  1  and  Remark  1  in  Renyi  (1967)  we  need  only  prove  that  for  any 
disjoint  intervals  Ej,  i  =  l,...n  1a(Ei)  i  =  !>—>n  are  independent.  If  this  is  proven,  it  then  follows 
from  the  argument  there  that  (N(t)  t  >  0}  is  a  process  with  independent  increments  and  consequently 
(N(t)  t  >  0}  is  Poisson  process. 

We  can  arrange  the  EjS  according  to  their  order  in  R+,  relabeling  them  if  necessary,  and  cut  them 
into  small  intervals  of  length  £  beginning  at  Ej.  For  any  £  >  0  there  exist  a  set  of  intervals 

Ej  i  =  1 . n  such  that  Ej  cEj,  )  Ei )  =  Cj£  Cj  are  integers  and  |E;|  -  |Ej  |  <  2e.  Furthermore, 

the  Ej  can  be  chosen  so  that  Fj  the  space  between  Ej  and  Ei+1  are  intervals  of  length  |  Fj  |  =  dj£, 
where  dj  are  also  integers,  i  =  l,...,n— 1  as  shown  in  the  diagram.  Denote  the  interval  of  length  £ 
which  contains  Fj_j  ^  Ej  by  Ait  and  the  interval  of  length  £  which  contains  Fj  f'j  E,  by  Ai2.  When 
Fj_j  f-'j  Ej  or  Fj  E;  is  empty,  define  A;;  =  empty.  Thus  |  Ay )  <  £,  j  =  1,2,  i  =  l,...,n. 


Diagram  1 


By  Lemma  1,  all  contiguous  subintervals  I  f  of  length  £  have  the  property  that  1a(|  are  indepen¬ 
dent.  It  follows  that  1a(eVs  are  independent  since  E1  Ej  =  empty  for  i  *  j  and  1A(E  >  only  depend 


on  1 


A(^)'s  such  that  Ij  cEr  Since  E,  c  E,  c  A,|  (j  E,  (j  A,:.  we 


have  for  any 


>.v  ,v:  /.v;  ,-.v 


<  ilf  .  .  .  ,im  <  n  m  = 

E  1  A(E{!)  UiEj 

*  E  ^  A(E1))  -  1  A(E1ib) 

m 

2:  E  fl  1  a(a,  ,  U  Ei,'  U  Ai,2) 
j=i 


=  E 


m  m 


,n  U(E-)  n  l1A(Ai  l)  1  A(A,  a)] 

b=>  j=> 


=  E- 

[  m  1 

n 1  a(e,  )  | 

[E- 

[a  ,  i 

ill  lA(A,,)  | 

[E- 

[n  '*<v] 

b=>  JJ 

i 

U=‘  J 

1 

b=I  J 

=  Esn  !a(e  •)  fe 

i=l  j 


-2Xm£ 


fhus 


'  m  1 

n 1  a(e,  )  | 

L  e'2Xm£ 

IH  ‘  J 

1 

m 

m 

n 1  a(e,  ) 

-nE 

j=l 

j=l 

m 

fj" 

"  fl  E  lA(E,) 


j=l 


- E  n 1  a(e,  ) "  -jn E  ^(Ei')  fe 
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-2XmE 


Since  £  are  arbitrary,  we  obtain  that  the  1a(EJ's  are  indePen^ent  from  ^  drt’  1  ndepen 

dent. 


For  the  nonhomogeneous  case,  we  have  similar  results. 


Theorem  T  Let  { N ( t )  t  >0}  be  a  counting  process  and  let  X  be  a  nonatomic  measure  on  B  the  class 
of  Borel  sets  on  R*  For  intervals  I  =  (a,b|  c  R\  let  X(I)  =  m(b)  -  m(a),  where  m  is  a  nonnegative 
nondesicasing  continuous  function  with  m(0)  -  0.  If  the  following  conditions  i)-iii)  are  satisfied,  then 
Niti  is  a  nonhomogeneous  Poisson  process  with  mean  function  E(N(t))  =  m(t). 

hi  For  an>  internal  I  =  ta.b),  P { N ( b )  -  N(a)  =  0}  =  e  /v(E 

in  P{  N(  h )  -  N ( a )  =  Fl }  >  2  }  =  o  ( ).( I )!  as  /.( I )  -4  0  uniformly  in  I. 


mi  t  or  anv 


timte  number  of  contiguous  intervals  I,  -  t,|  i=l,...,k  such  that 


/.il.  I  =  6j,  6;  >  0,  1=1 _ k,  the  rv  s  1  A(t ,  are  exchangeable. 


Proof  We  follow  the  same  line  of  argument  as  in  Lemma  1  and  Theorem  2.  We  need  only  replace 
p  -  e  ,'1"'  by  p  =  e  V.  divide  the  E,'s  by  sub  intervals  I,  with  X,(l,)  =  £|,  etc.  We  omit  the  details. 

Remark  I  Without  t)  and  u  >  in  Theorem  3,  the  condition  of  independent  increments  and  the 
exchangeability  of  tne  I  A,|  ,  in  ml  of  Theorem  3  are  not  equivalent. 

Remark  2  Definition  1,  the  condition  in  Theorem  1  (i.e.,  ^(E)  is  Poisson  distributed)  and  the  condi¬ 
tions  in  Theorem  3  are  equivalent.  We  show  this  in  the  following. 

Definition  1  implies  the  conditions  in  Theorem  3,  j.e.,  if  N(t)  has  independent  increments  then 


E  lAr,  c  l 


E  lA(i  c'  E  L 


All,  )  1  A(l,  )  -  c  ‘Ad,  )  ‘A(I,  ) 


=  e~m5'  =  E  1A„,-1 


A(I  i)  ‘  A(Im) 


for  any  1  <  i|  im  <  k  so  that  1  Aflj's  are  changeable. 


The  conditions  in  Theorem  1  imply  the  conditions  in  Theorem  3,  i.e.,  let  E  =  I, ,  E £  .  Then 

J=l 

^  1  At!,  i"  '  Ad,  ) 


=  P(^(I, )  =  0.  j  =  1 . m) 


P(^(E)  =  0)  =  e  >-<E)  =  e  mSl 


V  J  J  s 
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-  E  Ua.r'Ua.) 

for  any  1  <  ij  im  <  k.  Thus  the  l^jj's  are  exchangeable. 

By  Theorem  3  we  know  that  these  three  sets  of  conditions  are  equivalent.  However,  the  condi¬ 
tions  in  Theorem  3  are  somewhat  easier  to  check. 

Remark  3.  Oakes  (1972)  mentioned  the  result  that  if  £(1;),  i  =  where  the  I;  are  contiguous 

intervals,  are  independently  distributed  as  a  Poisson  distribution  with  mean  |I,|  for  all  m,  then  N(t) 
must  be  Poisson.  This  fact  also  follows  from  Theorem  2.  To  see  this  choose  special  contiguous  inter¬ 
vals  Ij  such  that  |Ij|  =8  then  £(Ij)  are  i,i,d  Poisson  (8).  Thus  the  1  a(I,)  s  are  exchangeable  for  any 
I-  i  =  of  length  8  and  contiguous.  By  Theorem  2  N(t)  is  a  Poisson  process. 

Remark  4.  Under  condition  iii)  of  Theoem  2  or  Theorem  3  we  can  weaken  conditions  i)  and  ii).  We 

need  only  assume  they  hold  for  the  special  class  of  intervals  having  the  form  I  =  (0  t], 

3.  Characterization  Via  Moments 

We  recall  some  definitions  and  properties  for  the  rth  joint  cumulant  which  was  discussed  in  a 
recent  paper  by  Block  and  Fang  (1985). 

Definition  2.  The  rth  joint  cumulant  of  (X  j,  .  .  .  ,  Xr)  is  defined  by 

”  A 

cum  (X„  •  ■  •  ,Xr)  =  £  (-lr1  (P-D!  Efl  Xj  -  En  Xj  (5) 

jevi  jevp 

c  J  J 

where  the  summation  extends  over  all  partitions  (Vj,  .  .  .  ,vp),  p  =  l,2,...,r,  of  { l,...,r}. 

We  state  some  properties  which  are  easy  to  check. 

(i)  If  any  group  of  the  X's  are  independent  of  the  remaining  X's  then  Cum  (Xj . Xr)  =  0. 

(ii)  For  (X[,  .  .  .  ,  Xr)  independent  of  (Yj,  .  .  .  ,Yr)  Cum  (X]  +  Yj,  .  .  .  ,  Xr  +  Yr) 

=  Cum  (Xj,  .  .  .  ,Xr)  +  Cum  (Y,,  .  .  .  ,Yr). 


(iii)  Cum  =  EXj  and  Cum  (Xj,  Xj)  =  Cov  (Xj,  Xj). 


(iv)  The  joint  cumulant  is  a  multilinear  operator. 


A  useful  relationship  between  moments  and  cumulants  is  given  in  the  following  lemma. 


Lemma  2.  If  E  |  Xj  |  m  <  °°  for  any  ; 


EXfXj,,  -  EXj— EXm  =  £  cum  (Xk,  kev,)—  cum  (Xk,kevp)  (6) 

where  ^  extends  over  all  panitions  (Vj,  .  .  .  ,Vp),  p  =  l,...,m-l  of  {l,...,m}. 

Proof.  See  Lemma  1  of  Block  and  Fang  (1985). 

For  the  characterization  problem  here  we  need  the  following  result. 

Lemma  3.  If  for  all  positive  integers  qj  <  rs  i  =  l,2,...m  such  that  £q,  <  J/j  we  have 
EYiq>  •  •  ■  Yj"  =  EY?1  ■  •  •  EYj-,  then  EY,r'....Y^  -  EY^  •  •  •  EY^m 


=  cum  (Ylt  .  .  .  ,  Y . . Ym,  .  .  .  ,  Ym) 

fi  u 


Proof.  By  Lemma  2. 


EY,r'  •••  Yi-  -  EY^  •  •  •  EY^ 


=  ^  cum  (Y,1,  iev,)-cum  (Y,\  ievp) 


where  the  sum  is  overall  partition  (Vp  .  .  .  ,vp)  p  =  l,2,...,m-l  of  { l,...,m}. 


Now  by  the  conditions  of  Lemma  3  and  induction  we  have  the  previous  sum  is 


cum  (Y,\  .  .  .  ,Y^) 


Using  Theorem  2.3.2.  of  Brillinger  (1975)  this  turns  out  to  be 


•  ■  t  *  *  «  »  i 


x  cum  (Yj,  jevjV-  cum  (Yj,  jevp)  (10) 

where  the  summation  is  over  all  indecomposable  partitions  of  the  set  wrt  to  parti¬ 
tion  where  the  number  of  i’s  is  r;  for  i  =  The  desired  result  follows 

by  using  the  conditions  of  Lemma  3  and  induction  once  again.  Since  for  all  proper  indecomposable 
partitions  Vj,  cum  (Yj,  j£V;)  =  0. 

Now  we  return  to  the  nonhomogeneous  Poisson  processes.  Some  basic  facts  are  reviewed  below 
for  the  nonhomogeneous  Poisson  process  (N(t),  t  >  0}  with  mean  m(t)  =  EN(t),  t  >  0. 

(i)  The  characteristic  function  of  N(t) 


<t>N(t)  (u)  =  E{exp  (iuN(t))}  =  exp  {m(t)(eiu-l)}  (11) 

is  analytic  in  the  whole  plane. 

(ii)  The  mth  moments  of  N(t)  are  polynominal  in  (J.  =  m(t).  More  precisely  we  have 


E{(N(t))k}  =  £  Si>k  ^ 

i=0 


(12) 


where  Sj  k  is  the  Stirling  numbers  of  the  second  kind  defined  by  formula 


%  =  -k  i  <-»H 

*•  j=0 


r 

i 


(13) 


(See  Riordan,  1937) 

We  characterize  the  nonhomogeneous  Poisson  process  as  follows: 

Theorem  4  Let  {N(t),  t  >  0}  be  a  counting  process  with  finite  moments  of  all  orders.  This  process  is 
a  nonhomogeneous  Poisson  process  if  and  only  if  for  all  real  tj,  i=l,...,k 

cum  (N(t,),...N(tk))  =  min  EN(tj)  .  (14) 

l  <  i  <  k 


Proof  If  (N(t),  t  >  0}  is  a  nonhomogeneous  Poisson  process,  we  assume  tj  <  t2  <  '<  tk.  Since  N(t) 
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has  independent  increments,  N^)  is  independent  of  N(tj)  -  N^),  i  =  2,...,k.  Thus  (N(t1),...,N(t1)) 
is  independent  of  (0,  N(t2)  -  N(ti),...,N(tk)  -  N(tj)).  Using  properties  of  the  cumulant 

cum  (NCtj)....^^)) 

=  cum  (N(tj),  N(tj)  +  NCt^-NCt,)....^^)  +  N(tk)-N(tj)) 

=  cum  (N(t1),...,N(t1)) 

+  cum(0,  N(t2)-N(t1),...,N(tk)-N(t1)) 

=  cum(N(t1),-.,N(t1))  (15) 

Let  C,k  denote  the  number  of  ways  of  partitioning  into  i  groups  for  where  i  >  k  or  i  <  1  we 

let  Ci>k  =  0.  For  example  when  k=2  we  have  C02  =  0>  C1>2  =  1,  C22  =  1,  C32  =  0  etc.  For  all  k, 
Ci  k  =  1,  Ckjk  =  1.  We  also  have  following  recurrence  formula. 

Q,k+1  =  i  Q,k  +  Q-l,k  (16) 

This  result  follows  since  one  can  count  the  number  of  ways  of  partitioning  { l,...,k+l }  into  i  groups  by 
considering  two  cases.  One  case  is  obtained  by  setting  one  element  aside,  then  partitioning  the  remain¬ 
ing  k  elements  into  i  groups  and  then  replacing  the  single  element  in  any  of  these  groups.  The  total 
number  is  iCi  k.  The  second  case  is  obtained  by  counting  the  numbers  of  partitioning  into  i  groups 
where  the  single  element  forms  one  group  and  then  partition  the  remaining  k  elements  into  i-1  groups. 
This  number  is  k. 

Notice  that  (15)  coincides  with  the  recurrence  formula  for  Sjk  in  (12)  and 

50, k  =  0  =  C0  k 

51, k  =  (-D1-1  [j  )  lk 

=  1  =  Cl,k 

So  we  get 

Si  k  =  Ci>k  .  (17) 

Now  using  this  fact,  we  can  prove  that  cum  (N(t1),...,N(t1 ))  =  jd,  where  |i  =  EN(t)).  We  use 
induction  on  k,  the  number  of  N(t;)'s  in  cum  (N(ti),...,N(t!)).  When  k=l  the  result  is  trivial.  For 


mm 


k=2,  Cov  (NCtj),  N^)  =  Var  N(tt)  =  JI,  which  is  well  known.  For  k  >  2,  by  Lemma  2  we  obtain 

cum 

=  E  {N(t!))k}  -  £  cum  (N(tj),  iev^-cum  (N(t,),  iev^  =  ti  (18) 

where  the  summation  ^  extends  over  all  partitions  V; . vp  of  for  p  —  2,...,k.  Now  by 

the  induction  assumption  and  (12),  (17) 

cum  (N(t1),..,N(t1)) 

=  E{(N(tl))k}  -  £  Ciik 

i=2 

=  icun1-icuni  =  cun  =  n  <19> 

i=l  i=2 

Combine  (15)  and  (19),  (14)  follows. 

To  prove  the  sufficiency  we  need  to  prove  for  any  interval  (tj,  t2],  N(t2)— N(ti)  is  Poisson  distri¬ 
buted  and  N(t)  has  independent  increments.  From  (14)  and  E[(N(t))k]  <  °°,  we  have 

cum  (N(t),...,N(t))  =  EN(t)  (20) 

Let  (J.  =  EN(t)  and  substitute  (20)  in  (18)  (replacing  tj  by  t),  we  have 

k  k 

E{N(t)k}  =  Vi  +  £  Ci,k  =  Z  c.,k  H* 
i=2  i=l 

which  coincides  with  the  moments  of  Poisson  rv’s  given  in  (12). 

It  is  well  known  (see  Y.  S.  Chow  1978,  p.  280)  that  if  the  characteristic  function  of  a  distribution 
F  is  analytic,  the  F  is  uniquely  determined  by  its’  moment  sequence.  Since  the  characteristic  function 
of  Poisson  is  analytic  in  the  whole  plane,  we  have  that  N(t)  is  Poisson  with  mean  E(N(t)). 

For  increments  N(t2)  -  N(t[)  on  (t^  t2]  we  need  only  check 

cum  (N(tj)  -  N(tj),...,N(t2)  -  N(tj))  =  E(N(t2)  -  N(t,))  (21) 

since  the  above  argument  gives  that  N(t2)  —  N(t[)  is  Poisson.  By  (14)  it  is  easy  to  see  that 
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cum  (N(tj)  -  NCtj),...^)  -  N(tj)) 

=  cum(N(t2),  N(t2)  -  N(t1),...fN(T2)  -  N(t,)) 

-  cum(N(t1),  N(t2)  -  N(t,),...,N(t2)  -  N(tO) 

=  cum(N(t2),  N(t2)-N(t0,-Mt2)  -  N(ti» 

=  cum(N(t2),...,N(t2),  N(t2)  -  NdO) 

=  EN(t2)  -  EN(tj)  =  E(N(t2)  -  N(t,)) 

Here  we  use  the  multilinear  property  of  joint  cumulant. 

It  follows  from  (21)  that  N(t2)  -  N(tj)  is  Poisson  distributed  with  mean  E{N(t2)  -  N(t j ) }.  Con¬ 
sequently,  the  characteristic  function  of  N(tj)  -  N(ti_1)  is  analytic.  We  also  know  that,  for  rv’s 
Xj,  .  .  .  ,Xr  if  for  each  i,  i  =  l,...m 

I  EX,1  (22) 

has  a  positive  radius  of  convergence  as  a  power  series  in  8,  then  the  joint  distribution  of 

(X1?  .  .  .  ,Xn)  is  determined  by  its  cross  moment  EXj1,  .  .  .  ,  X ^  where  r;  are  nonnegative  integers 
(see  Billingsley,  1979,  p.  351)  for  Poisson  marginal  (22)  hold.  To  prove  independent  increments  we 
only  need  to  check 

EX1\...fXi-  =  EX[*,...,EXi-  (23) 

for  all  Tj,  i=l,...m  where  Xj  =  N(tj)  -  N(tj_j)  and  where  (tj_j  tj,  i  =  are  disjoint  intervals  on 


We  check  (23)  by  using  Lemma  3  and  induction.  For  m=l,  its’  trivial.  For  m=2,  fj  =  r2  =  1, 

Cov  (Xlf  X2)  =  Cov  (N(tj)  -  N(to),  N(t, )  -  N(t,)) 

=  Cov  (N(tj)  N(t2))  -  Cov  (N(to),  N(t2)) 

-  Cov  (N(t,),  N(t2))  +  Cov  (N(t0),  N(t,)) 

=  EN(tj)  -  EN(to)  -  EN(t,)  +  EN(to)  =  0 

So  EXi  X2  =  EXj  EX2.  Similarly,  we  have 

Cov  (X,  Xj)  =  0  i  *  j 


WWW! 
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which  means  increments  on  disjoint  intervals  are  uncorrelated 
By  (14),  we  also  have 

cum  (N(t,)  -  N(t0),...,N(tm)  -  NCV-,))  =  0  (25) 

which  means  increments  on  disjoint  intervals  have  zero  joint  cumulant.  More  generally, 

-Jk  Jm 

cum  (X, . X, . X,; . XJ 

=  2r~'  (EN(tj)  -  EN(tj)} 


+ 


I 

2H 


2r®  i 


{EN(to)-EN(to)}=0 


(26) 


We  do  a  simple  case  as  example.  Assume  <  t(  <  t2. 


cum  (Xj,  Xj,  X2) 

.=  cum  (N(ti)  -  N(to),  N(tj)  -  N(to),  (N(t2)  -  N(t,)) 

=  cum  (N(t,),  N(t,),  N(t2)  -  cum  (N(t,),  N(t,),  N(t,)) 

-  cum  (N(t,),  (N(to),  N(t2))  +  cum  (N(t,),  N(to),  N(t,)) 

-  cum  (N(to),  N(t,),  N(t2))  +  cum  (N(to),  N(t,),  N(t, )) 

+  cum  (N(to),  N(to),  N(t2)  -  cum  (N(t0),  N(t0),  (N(t, )) 

=  { EN(t, )  -  EN(t, ) }  +  3  {EN(to)  -  EN(t0)}  =  0 

From  (26),  Lemma  3  and  induction  (23)  follows  for  all  r,,  which  imply  independent  increments  for 


N(t).  This  completes  the  proof. 
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As  a  by-product  we  have  the  following  result. 


Theorem  5.  If  the  rv  X  satisfies  E  |  X  | n  <  <»  for  all  n,  the  X  is  Poisson  distributed  if  and  only  if 
cum  (X,...,X)  =  EX  for  all  m. 


The  characterization  here  provides  another  example  of  independence  characterized  via  uncorrela¬ 


tion.  Test  methods  for  a  counting  process  being  Poisson  can  be  constructed  by  using  this  result. 
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